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Abstract 

In light of the multiple sex chromosome systems observed in howler monkeys (Alouatta Lacepede, 1 799) 
a combined cladistic analysis using chromosomal and molecular characters was applied to discuss the 
possible origin of these systems. Mesoamerican and South American howlers were karyologically com- 
pared. FISH analysis using the chromosome painting probes for the #3 and #15 human chromosomes 
was applied to corroborate the homeology of the sexual systems. We found that the HSA3/15 syntenic 
association, present in the sex chromosome systems of South American Howlers, is not present in those 
of Mesoamerican ones. The autosomes involved in the translocation that formed the sexual systems in 
the Mesoamerican and South American species are different, thus suggesting an independent origin. 
Parsimony analysis resolved the phylogenetic relationships among howler species, demonstrating utility 
of the combined approach. A hypothesis for the origin of the multiple sex chromosome systems for the 
genus is proposed. 
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Introduction 

Howler monkeys (genus Alouatta Lacepede, 1799 of the family Atelidae) exhibit one 
of the widest geographic distributions recorded to date for Neotropical Primates. 
Their distribution extends from southern Mexico to northern Argentina (Crockett 
and Eisenberg 1987, Rylands 2000). They inhabit a diverse range of environments, 
including tropical rain forests, flood forests, gallery forests, patches of forest and de- 
ciduous and semideciduous seasonal environments (Crockett and Eisenberg 1987, 
Zunino et al. 2001). There remains a lack of consensus regarding both the number 
of species within the genus, which, depending on the author, ranges from 9 to 14 
species (Rylands 2000, Groves 2001, 2005, Gregorin 2006, Rylands and Mittermei- 
er 2009), and the phylogenetic relationships among them. This shows the complex- 
ity of the taxonomy of Alouatta and highlights the importance of including a larger 
number of variables for a more accurate characterization of the species in the genus. 
We adhere to the classification proposed by Groves (2001, 2005) in recognizing 10 
species {Alouatta belzebul Linnaeus, 1766, A. seniculus Linnaeus, 1766, A. sara El- 
liot, 1910, A. macconnelli Linnaeus, 1766, A. caraya Humboldt, 1812, A. palliata 
Gray, 1849, A. pigra Lawrence, 1933, A. guariba Humboldt, 1812, A. nigerrima 
Lonnberg, 1941, A. coibensis Thomas, 1912), since it considers both morphological 
and genetic information. 

In this context, and to contribute to the description of the phylogenetic relation- 
ships in the genus, several authors have proposed that chromosomal data can also be 
used as phylogenetic markers, since they are inherited as mendelian characters and are 
conserved within species (Sankoff 2003, Dobigny et al. 2004, Stanyon et al. 2008). 
Following the Maximum Parsimony criterion, karyological comparisons allow the 
identification of chromosomal forms shared by common ancestrality. 

In primates, different researchers in the last three decades have proposed chro- 
mosomal speciation as a probable evolutionary mechanism to explain the diversity 
observed in living species (de Grouchy et al. 1972, Seuanez 1979, Dutrillaux and 
Couturier 1981, Clemente et al. 1990, Stanyon et al. 2008, de Oliveira et al. 2012). 
In howler monkeys, species exhibit diploid numbers (2N) ranging from 44 in Al- 
ouatta seniculus to 58 in Alouatta pigra, and in a large number of species, multiple sex 
chromosome systems in males originated from Y-autosome translocations have been 
described (Table 1). The chromosomes involved in the Y-autosome translocations in 
Alouatta caraya, Alouatta macconnelli, Alouatta guariba guariba Humboldt, 1812, Al- 
ouatta guariba clamitans Cabrera, 1940, Alouatta sara and Alouatta seniculus arctoidea 
Cabrera, 1940, are homeologous to the same regions of human chromosomes #3 and 
#15 (Consigliere et al. 1996, 1998, Mudry et al. 2001, de Oliveira et al. 2002). 

The phylogenies proposed so far for Alouatta have used either molecular markers 
(y'-globin (Meireles et al. 1999), Mt ATP synt 8 and 6, Mt cyt b, CAL and RAG1 
(Cortes-Ortiz et al. 2003)) or chromosomal characters (de Oliveira et al. 2002). 
However, the combination of different variables can improve the phylogenetic sig- 
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Table I. Cytogenetic characteristics of howler monkeys {Alouatta). 



Species 


2N 


Sex Chromosome Systems 


References 


A. belzebul 


950 c?49 


X 1 X 1 X 2 X 2 /X 1 X 2 Y 


Armada et al. 1987 § 


A. s. seniculus 


9c?47 to 49t 


XY 


Yunis et al. 1976 


A. s. stramineus 


?c?47 to 49t 


X,X,X J X,/X,X ; Y,Y, 


Lima and Seuanez 1991 5 


A. s. arctoidea 


$44 c?45* 


X,X,X r !C/X,X,Y,Y, 


Stanyon et al. 1995 


A. sara 


$c? 48 to 51t 


X.X.WW 


Minezawa et al. 1985 




9c?50 


XXX X / XXY.Y, 

112 2 1^12 


Stanyon et al. 1995 


A. macconnelli 


9 c? 47 to 49t 


X,XXX/X,XYY 

11^2 1212 


Limaetal. 1990 


A. caraya 


9<?52 


XX/XY 


Egozcue and De Egozcue 1966, 
Mudryetal. 1984, 1994 






X,X,X^,/X,X,Y,Y, 


Rahn et al. 1996 s , Mudry et al. 1998 5 , 2001 5 


A. palliata 


9c?56 


XX/XY 


Torres and Ramirez 2003 




954 c?53 


X 1 X 1 X 2 X 2 /X 1 X J Y 


Maetal. 1976 
Solari and Rahn 2005 s 


A. pigra 


9c?58 


X,X,X 2 X 2 /X 1 X 2 Y 1 Y 2 


Steinberg et al. 2008 s 


A. guariba guariba 


950 c?49 


XX/XY 


Koiffmann and Saldanha 1974 




c?49 


X,X,Y 


de Oliveira et al. 1995 




950cj49 


X | X 1 X,X,X,X,/X,X,X,Y,Y, 


de Oliveira et al. 2002 


A. guariba clamitans 


946 c?45 


XX/XY 


de Oliveira et al. 1995 






X l X l X 2 X 2 /X l X 2 Y 


de Oliveira et al. 1998 






x,x,x,x,x,x,/x,w,y, 


de Oliveira et al. 2002 


A. nigerrima 


950 


XX 


Armada et al. 1987 


A. coibensis 


ND 


ND 





tThese differences are due to the presence of microchromosomes (1 to 3 per nuclei); % Differences due 
to a variation in microchromosome number between sexes. ND: not yet cytogenetically characterized. 
§ Meiotic studies performed to corroborate the sex chromosome system. 



nal due to the possible common shared history of different datasets. At the same 
time, this combination can increase the support of a tree, since different characters 
evolve at particular rates and will support different parts of the tree (Kluge 1989, 
Whittaker et al. 2007). 

In the present contribution, howler species were karyologically compared and 
FISH analyses were carried out to corroborate the homeology of the sex chromosome 
systems among them. Using these data and molecular data obtained from the litera- 
ture, a phylogenetic analysis combining them in a single matrix was performed. 



Methods 

Sampled specimens: A total of 29 adult specimens of both sexes of four species of howl- 
ers, both from captivity as well as from the wild within their natural geographical 
distribution, were analyzed: Alouatta caraya (9 S and 6 A. guariba clamitans (1 S), 
A. pigra (6 $ and 5 $) and A. palliata (2 S). 
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The origin of the animals was as follows: 
Argentina 

A. caraya, 1 S from Corrientes Zoo, Corrientes; 1 S from Ecological Park "El Puma", 
Misiones; 1 o and 2 9 from Mendoza Zoo, Mendoza; 6 o and 4 9 from the Black 
Howler Monkey Reeducational Center, La Cumbre, Cordoba. 
A. g. clamitans, 1 S from "Giiira-Oga", Misiones. 

Mexico 

A. pigra: 4 $ and 4 $ were sampled in the wild in Campeche, Yucatan Peninsula; 2 S 

and 1 $ from San Juan de Aragon Zoo, Mexico City. 
A. palliata: 1 S from San Juan de Aragon Zoo, Mexico City; 1 S from Chapultepec 

Zoo, Mexico City. 



Classical cytogenetic analysis 

Chromosome preparation: Peripheral blood samples were collected from all animals with 
previously heparinized disposable syringes. Lymphocytes were cultured for 72 h at 37 
°C following Mudry (1990). At least 50 metaphases were analyzed to determine the 
diploid number (2N) at 1000x. Metaphase spreads were treated with G-Wright band- 
ing (Steinberg et al. 2007). At least 10 G-banded metaphases with the species diploid 
number (2N) were photographed with a Leica DFC 340 FX camera. Chromosomes 
were arranged according to previously described karyotypes using Photoshop CS (Ado- 
be) and the species assignation of each specimen was corroborated. 

Analysis of homeologies: For A. caraya and A. g clamitans, the homeologies with 
human chromosomes and the homeologies with the other South American howlers 
are well known (Consigliere et al. 1998, Mudry et al. 2001, de Oliveira et al. 2002, 
Stanyon et al. 201 1). The G-banded chromosomes oi A. pigra and A. palliata were first 
compared with those of A. caraya and A. g. clamitans. We took A. caraya % karyotype 
as the reference for the comparisons with Mesoamerican howlers (Mudry et al. 200 1 , 
Szapkievich and Mudry 2003). To compare homeologies, the G-banded metaphases 
obtained for A. caraya, A. g. clamitans, A. pigra and A. palliata were also compared with 
those published for A. g. guariba (de Oliveira et al. 2002, Stanyon et al. 201 1), A. mac- 
connelli (de Oliveira et al. 2002), A. s. arctoidea (Consigliere et al. 1996), A. belzebul 
(Armada et al. 1987, Consigliere et al. 1998) and A. sara (Consigliere et al. 1996). 



Cytomolecular study 

FISH analysis with human chromosome painting probes #3 and #15 was used as a tool 
to confirm the identity of the sex chromosome systems in howlers. Whole chromo- 
some painting probes for human chromosomes #3 (red), #15 (green), #21 (green), X 
(green) and Y (red) (PCT3 Cy3, PCT15 FITC, PCT21 FITC, PCTX FITC, PCTY 
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Cy3, LEXEL S.R.L., Buenos Aires, Argentina) were used for FISH analysis on the 
metaphases of A. pigra, A. caraya, A. g. clamitans and A. palliata. Homo sapiens (HSA) 
metaphases were used as a positive control of hybridization. The HSA3/21 syntenic 
association, considered ancestral in mammals and conserved in most primate species 
(Miiller et al. 2000), was analyzed simultaneously as a control of synteny conservation. 
Human X and Y chromosomes were also tested. 

FISH was performed according to the supplier's instructions (LEXEL S.R.L., Bue- 
nos Aires, Argentina). Slides were counterstained with DAPI (Sigma) and analyzed 
with a Leica DMLB fluorescence microscope. Chromosome images were obtained 
with a Leica DFC 340 FX camera. Images were processed with Image Pro-Plus 4.5 
(Media Cybernetics Inc.). 

Our results were compared with those previously described (Consigliere et al. 
1996, 1998, Mudry et al. 2001, de Oliveira et al. 2002, Stanyon et al. 201 1). 



Phylogenetic analysis 

Chromosomal dataset: We used data obtained from the comparisons of G-banding pat- 
terns and the analysis of chromosomal syntenic associations, both from the present 
study and from previous reports (Consigliere et al. 1996, 1998, Garcia et al. 2001, 
2002, Mudry et al. 2001, de Oliveira et al. 2002, Amaral et al. 2008, Stanyon et al. 
2001, 201 1). We considered the structural changes as characters. The pattern observed 
before and after their occurrence, i.e. their presence or absence, was considered as the 
character states. The matrix (see Appendix 1) was produced taking into considera- 
tion the characters proposed by Neusser et al. (2001) and modified for howlers by de 
Oliveira et al. (2002). These authors used an abbreviated nomenclature for ancestral 
Platyrrhini chromosome forms with a correspondence in the human karyotype. In the 
present contribution, for the character nomenclature, we refer directly to the human 
G band ideogram (Table 2). New characters were obtained from our karyological com- 
parisons and introduced in the chromosomal dataset. 

Molecular dataset: The sequences available in GenBank for the same species used in 
the G-banding pattern and FISH comparisons were taken into to choose the molecular 
marker. The only molecular marker that fullfiled all the requirements was cyt b. The se- 
quences used were (Genbank Accession Numbers): A. helzehul (AY374348.2), A. caraya 
(AY374359.2), A. s. arctoidea (AY065886.1), A. sara (AY065887.1), A. macconnelli 
(AY065888.1), A. g. guariba (AY065899.1), A. g. clamitans (DQ679782.1), A. pigra 
(AY065884. 1), A. palliata (AY065879. 1) (Cortes-Ortiz et al. 2003, Harris et al. 2005, 
Lorenz et al. 2005, Nascimento et al. 2005, Casado et al. 2010). Cebus apella Linnaeus, 
1758 (FJ529102.1) and Lagothrix lagotricha Humboldt, 1812 (AY671799.1) were used 
as outgroups. C. apella, from the Cebidae family, was taken as an outgroup species, since 
it is accepted that this species presents the most ancestral karyotype within Platyrrhini 
(Clemente et al. 1990, Garcia et al. 2000). Lagothrix lagotricha, also a member of the 
Atelidae family, was chosen as the second outgroup to test the monophyly of the group. 
All sequences were aligned using CLUSTALW (Thompson et al. 1994). 
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Table 2. Human chromosome syntenic association considered as characters and used to construct the 
binary matrix of chromosomal homeologies among howler monkeys (modified from Neusser et al. 2001, 
de Oliveira et al. 2002). 



1 


1 n91 nrpr/1 nl 9 91 

i uz i -pici/ 1 p i z-z i 


3 A 


1 Snl 1 nl 3- n9S nrer/V 
ijqii ui j) qz j qici/ l 


Uy. 


(\ On/1 AlV) 

U u q/ lt} pJ 2 


2 


Sn^l H nrpr/ 7n99- n! 1 n91 

jqj i .j-qiei/ /pzz, qii qzi 


37 
j / ■ 


Vnfpr n 1 9 
ZpLCl U1Z 


70 


I 1 H\ On/1 Artl 

I I / ^iuq/ iop^ 7 


2 

J. 


5pter-q3 1 .2/5q3 1 .3-qter 


JO. 


1 An 

i oq 


71 

/ 1 . 


1 On/1 An 

iuq/ lop 


4 


9 n r nl 7/1 hn 

zpucr qiZ/ lOq 


.7./. 


3n94 n9 1 ■ n 1 3 n9A/ 1 Sn 1 1 

jpzi qz i , q i j qzo/ i jq 1 1 


79 
/ z. 


1 On/1 (in/4ni-pr n99 
iuq/ lop/ipuci-qzz 


j ■ 


4n^l 3 nrpr/4n93 n^l 9 
iqj 1 . J-qLCl/ iqZ J-qj 1 .Z 




ql3; q25 qter 


7-t 
/ J. 


1 Ort/1 On/1 An/4nrpr n99 

iup/ iuq/ lop/ipcer-qzz 




4n93 n^l llArttpr n99 

iqz j-qj i .z/ 'ipici-qzz 


4(1 


1 1 /Snrpr n^l 9 

11/ jptcr-qj i .z 


74 


1 On/1 On/1 An 

iup/ iuq/ lop 


7 


(^(^n/^ 


41 


S nf p r n^l 9/7 n 99- nl 1 ■ n91 

jpcer-qj i .z/ / pzz, qi i , qz i 


/ J- 


1 0/1 3 

1 J/ u 


g 


A 


42 


1 9/0 

1 Z/ y 


/ o. 


1 0/99 
i yi zz 


q 

7. 


Rn/1 £ 
op/ 1 o 


4} 


lp21-pter/2pter ql2 


77 
/ / • 


99/1 n91 ntpr 

zz/ ipzi-ptcr 


1 0 
1 U. 


1Sn91 3 n94/1Sn13 n91 9 
1 jqZ I.J-IJZI/ 1 jqi J-qZ 1 .z 


44 


i oq/ ipccr-qzz 


7S 
/ o. 


1 0/99/1 n9 1 nrpr 

i y i zz/ 1 pz i 'pier 


] I 


l^nll nl V n9S n t~i^ r 

Ijqil qij, C|Zj qtci 


4S 


22/14 


70 
/ J. 


9nl ^ nrpi74n93 n^l 9 

zqi j qicr/iqz j-qj i .z 


12 


7n91 nl 1 • nl 1 n91 ■ n99 ntpr 
/ pz 1 p 1 1 , (J 1 1 (JZ i , (JZZ CJCCI 


46. 


9nl 3 nrer/90 
ZqiJ (JLCl/ZU 


80 


A/1 nl 9 91 
o/ 1 p 1 Z-Z 1 


1 3 


sq 


47 


9n1 3 nrpr/4n9^-n31 9 
ZU1J (JLCi/ 4:UZJ-(J J 1 .Z 


81 


Y/1 Snl 1 nl ^- n9S nrpr/ ^n94 

1/ ijqii qu, qz j qrer/ jpzf 


1 4 


1 z 


48 


Xn/9nl ntpr 

oq/ zq i j qicr 




n9 1 • nl ^ n9A 

qz i , qi j qzo 


1 j. 


1 1 

1 i 


49 


7n99- nl 1 ■ n91 /8n 
/ pzz, 111 1, CJZ1/OC1 


H9 


1 nl 9 91 /Snrpi- n^l 9/7n99- 

ip iz-z 1/ jpcci-qj i .z/ / pzz, 


1 


1 D 


ju. 


7n99- nl 1 ■ n91 /8n 

/ pzz, q 1 1 , cjzi/oti 




nl 1 ■ n91 /Sn^l ^ nrpr/ 7n99- 
qi i , qz i / jqj i . j-qcer/ / pzz, 


1 7 


C) 

7 


S 1 
J i . 


17/9ntRr nl 9/1 9 

i / / zpccr q i z/ 1 z 




nl 1 n91 

q 1 1 qz i 


18 


*s n t-fi i-n94*n/l nl 7* /i 1 / nl V 

jpccr pzi, pz i piz, (J1Z (11 J, 


JZ. 


9r»fpr nl 9/1 9 

zpcer qiz/ iz 


Q 5 
o J. 


0/99 

J/ZZ 




q27 qter 


j J ■ 


1 n^9 nrpr/1 1 

i qjz qtci/ 1 1 


84 


1 7/1 1 

1//H 


1 0 


3n94 n91 ■ nl 3 n9A 

jpzi qzi, qi j qZO 


^4 


ln^9 nfpr/Cl 1 /^.nl-^r n^l 9^ 

iqjz qter/i^i 1/ jpter-qji.zj 2 


OJ. 


Jr-rpr n94- n91 nl9-nl9nl^' 

jprer pzi, pzi piz, qiz qi j, 


zu. 


lq32 qter 


j j ■ 


H 1 /Snrpr n^l 9^ 

l, 1 1/ jpccr-qj i .z^ 2 




n 9 7 n t~f* r I X n 

qz / qici / op 


9 1 


1 n91 n^l 

1 qz 1 q 9 1 


SA 

JO. 


1 R/l 4 
lo/ 11 


R£ 
ou. 


1 ^n9 1 ^ n94/l n^9 nrpr 

i jqz i . j-qzi/ iqjz qrer 


22 


^nl 9/91 
jp 1Z/Z1 


S7 

j / . 


Writer n94- r>9 1 r*19*n19n!3- 

jpccr pz^±, pzi piz, qiz qi j, 


R7 
o / . 


!An/lSn91 H n94/ln^9 nrpr 

ioq/ 1 jqz i . j-qzi/ iqjz qtci 


9 3 


1 (In 

lup 




n97 nter/1 Sn91 3-n94 
qz/ uLci/ 1 juz i . j'lizt 


88 


4nrpr n99/lnl9 91 

iprci-uzz/ 1 p i z-z i 


24 


22 


Jo. 


Hnrer r>94- n? 1 nl9-n19n!3- 

jptCl pZT:, pz 1 piZ, U1Z (Jl J, 


RO 
0 J. 


14 


9S 

z j. 


v 




n97 ntpr/1 Sn91 ^n94/l/^n 

qz/ qtcr/ 1 jqzi.j-qzi/ loq 


Oft 
JU. 


2pter ql2/4pter-q22 


7 A 
ZD. 


v 
I 


SO 
J J . 


1 Sn91 ^ n94/l An 
i jqz i . j-qz4:/ loq 


0 1 
7 1 . 


!Snl^n91 9/7n99-nll- 

i jqi j-qz i.x// pzz, qi i , 


27, 


Snrpr-nHI 9 


60. 


17/10p 




q21/5q31.3-qter/ 7p22; qll; 


28. 


(5q31.3-qter/7p22; qll 


61. 


17/10p/19 




q21 




q21) 2 


62. 


10p/19 


92. 


6/15q21.3-q24 


29. 


lp21-pter 


63. 


22/20 


93. 


14/Ipl2-21 


30. 


lpl2-21 


64. 


22/20/lq21 q31 


94. 


6/15q21.3-q24/l4/lpl2-21 


31. 


4q31.3-qter 


65. 


20/lq21 q31 


95. 


17/8p/18 


32. 


4q23-q31.2/15ql3-q212 


66. 


(ll/5pter-q31.2)^ 


96. 


22/10q/16p 


33. 


4pter-q22 


67. 


lpl2-21/8p 


97. 


2ql3qter/ll 


34. 


l4/15q21.3-q24 


68. 


7p21 pll; qll q21; q22 qter 


98. 


(10q/16p) 2 /lq21 q31/20 


35. 


(lOq/^pydOq/^p), 




/!4/15q21.3-q24 


99. 


Y/7 



/: separates the chromosomal segments that constitute an association. ( ) : n= number of repeats in the 
segment 



Phylogeny: A Maximum Parsimony phylogeny using the exhaustive search option was 
obtained with PAUP 4.0 software (Phylogenetic Analysis Using Maximum Parsimony, 
(Swofford 2002)), for each separate partition and the combination of both the chromo- 
somal and molecular datasets. All characters had the same weight, based on the premise that 
chromosome rearrangements occur by equal chance (de Oliveira et al. 2002, Dobigny et al. 
2004). The relative stability of nodes was assessed by bootstrap estimates (Felsenstein 1985) 
based on 200 iterations. Each bootstrap replicate involved a heuristic parsimony search 
with 10 random taxon additions and tree-bisection reconnection (TBR) branch swapping. 
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Results 

Classical cytogenetic analysis 

Karyological analysis: The cytogenetic characterization of the Alouatta specimens 
showed diploid numbers, sex chromosome systems and G-bandings patterns in agree- 
ment with the ones previously described for each species. Figures la and lb show all 
the comparisons performed. 

Chromosomal homeologies between A. caraya and A. palliata: The chromosomal rear- 
rangements that could explain the homeologies were grouped in two categories: 1) A. 
palliata chromosomes with no rearrangements with respect to A. caraya chromosomes: 
3, 5, 6, 7, 8, 9, 10, 13, 14, 16, 17, 19, 20, 21, 22, 24, 25 and X,; 2) A. palliata chro- 
mosomes with more than one rearrangement with respect to A. caraya chromosomes: 
1, 2, 4, 11, 12, 15, 18 and 23. No homeologies were allocated for A. palliata chromo- 
some 26 and chromosome arms 4q and 2p using the level of resolution of the classical 
cytogenetic techniques applied. The rearrangements detected between the A. caraya 
and A. palliata karyotypes included at least seven fissions/fusions, two paracentric in- 
versions and one deletion. A. caraya chromosome 7 (X, in males) shares homeology 
with two A. palliata chromosome pairs, 23 and 18, which are not the ones involved in 
the sex chromosome system in A. palliata. The A. palliata chromosomal pair 19 (3^ in 
males) shares homeology with chromosome 14 of A. caraya. 

Chromosomal homeologies between A. caraya and A. pigra: The chromosomal rear- 
rangements that could explain the homeologies were grouped in two categories: 1) 
A. pigra chromosomes with no rearrangements with respect to A. caraya chromo- 
somes: 2, 5, 6, 7, 8, 10, 15, 16, 17, 20, 22, 23, 25, 28 and X,; 2) A. pigra chromo- 
somes with more than one rearrangement with respect to A. caraya chromosomes: 
1, 3, 4, 9, 11, 13, 18, 19, 24, 26 and 27. No homeologies were allocated for A. 
pigra chromosomes 4p prox , 12, 14 and 21 using the level of resolution of the classical 
cytogenetic techniques applied. The rearrangements detected between the A. caraya 
and A. pigra karyotypes included at least 12 fissions/fusions, two paracentric inver- 
sions, two translocations and one deletion. A. caraya chromosome 7 (X 2 in males) 
shares homeology with two A. pigra chromosome pairs, 26 and 19, which are not 
the ones involved in the sex chromosome system in A. pigra. A. pigra chromosome 
17 (X 2 in males) shares homeology with chromosome 14 of A. caraya (which in turn 
has homeology with HSA7). 

Chromosomal homeologies among all howlers: The chromosomal homeologies found 
among all howlers are shown in Table 3 and Figures la and lb. Results show that 
Mesoamerican howlers share several human chromosomal syntenic associations with 
South American ones: HSA15ql3-q21.2/4q23-q31.2 and HSAl6p/10q, shared with 
all howlers; HSA15q21.3-q24/l4, shared with all howlers except A s. arctoidea and 
A. macconnelli, and HSA8p/18, shared with all howlers except A. s. arctoidea. Two 
new chromosomal syntenic associations, HSA4p ter -q22/9/l 1 and HSA15q21.3- 
q24/l4/21q, were found for A. pigra in chromosomes 1 and 4q, respectively. 
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Figure I . Comparison of Homo sapiens (HSA) , Alouatta caraya (ACA), A pigra (API), A palliata (APA) and 
A. guariba clamitans (AGUc) G-banded chromosomes, taking A caraya 's karyotype as reference. On the left, 
human chromosomal bands with homeology for its corresponding ACA chromosome segment are indicated. 
The boxes highlight the homeologies of the autosomes involved in the sex chromosome systems in these 
species a Comparison for ACA chromosomes #1 to #13 b Comparison for ACA chromosomes #14 to X, 
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Table 3. Chromosomal homeologies between howlers, obtained from data both from this contribution 
and from previous reports. ACA: A. camya; API: A. pigra; APA: A. palliata; AGU: A. guariba; ASEa: A. s. 
arctoidea; AMA: A. macconnelli; ASA: A. sara; ABE: A. belzebul. 
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Cytomolecular study 

In the Homo sapiens metaphases, the hybridization signals on chromosomes HSA3, 
HSA21, HSA15, HSAX and HSAY for chromosome painting probes #3 (red), #21 
(green), #15 (green) (Figure 2a), X and Y (data not shown) were corroborated. 

In A. g. clamitans, the signal for HSA21 was observed in 18q tcr , the signal forHSA3 
was observed in 18q prox (thus corroborating the HSA3/21 synteny in A. g. clamitans), 
I6q, 7q (X 2 in males) and Y , and the signal for HSA15 was observed in 1. , 2. , 7p 
(X^) and Y . This corroborates the HSA3/15 syntenic association to the multiple sex 
chromosome system X J X 1 X 2 X 2 X 3 X 3 /X 1 X 7 X 3 Y 1 Y 2 of this species (Figures 2b and 3b). 

In A. pigra, the signal for HSA3 was observed in 16q and 19q, while that for 
HSA21 hybridized in 4p , thus indicating that the HSA3/21 synteny is not present in 
A. pigra (Figures 2c and 3c). The probe for HSA15 hybridized in A. pigra metaphases 
in 4q Kr , 22q t r and 26q, showing that the HSA3/15 syntenic association is also absent. 
None of these A. pigra chromosomes is involved in the sex chromosome system of this 
species (Figures 2d and 3c). 

In A caraya, the signal for HSA21 was observed in 21q Kr , whereas that for HSA3 
was observed in 21q , thus confirming the conservation of the HSA3/21 synteny. 
HSA15 hybridized in 7p (X 2 in males) and Y lKr , and HSA3 in 7q and Y Iprox , exhibiting 
the HSA3/15 syntenic association in the sex chromosome system X XjXJC/X XYjY 
(Figure 3a). 

A. palliata showed a pattern similar to that of A. pigra (therefore Figure 2 illustrates 
only the latter). HSA3 hybridized in 16q and 18q, HSA21 hybridized in lp and 
HSA15 in 2q Kr , 21q ter and 23q (Figure 3d). Both the HSA3/21 and HSA3/15 syntenic 
associations are absent in A. palliata and chromosomes with homeology to HSA3 and 
HSA15 are also not involved in the sex chromosome system of this species. 

The probe for the human X chromosome showed positive hybridization signal in 
Xj of all the species analyzed. The probe for the human Y chromosome did not hybrid- 
ize in any of the howler species (data not shown). 
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Figure 2. Analysis of the conservation of the HSA3/21 and HSA3/15 syntenic chromosomal associations 
in howlers (bar=10 (im). The arrows indicate the chromosomes with positive FISH signal a Homo sapiens 
partial metaphase hybridized with probes HSA21 (green) and HSA3 (red) (control of the hybridization). 
Inset: Homo sapiens partial metaphase hybridized with HSA15 (green) b AGUc metaphase hybridized 
with HSA15 (green) and HSA3 (red) c API metaphase hybridized with HSA3 (red) and HSA21 (green) 
d API metaphase hybridized with HSA15 (green). 



Phylogenetic analysis 

The data obtained from the G-banding pattern and FISH homeologies, together with 
cyt b sequences obtained from previous reports, were used as the basis to perform a 
cladistic analysis. The HSAY/7 association, corresponding to the Y- autosome translo- 
cation that gave rise to the multivalents observed in A. pigra and A palliata, was added 
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d. 1 2 16 18 21 23 

Figure 3. Howler monkeys G-banded chromosomes with positive signal for the human chromosome 
painting probes analyzed. On the right, the hybridization pattern of human chromosomes #3, #21 and 
#15. a A. caraya b A. guariba clamitans c A. pigra d A. palliata. 



as an extra character to the original list (de Oliveira et al. 2002). The syntenic associa- 
tions HSA4p rer -q22/9/l 1 observed in chromosome 1 of A. pigra and HSA15q21.3- 
q24/l4/21q observed in chromosome arm 4q were not included in the analysis, be- 
cause, as autopomorphies for A. pigra, they are considered non-informative. 

Three data matrices were obtained: one including only chromosomal data, another 
including only molecular data and the last one including both types of characters 
(chromosomal and molecular) in a single matrix (see Appendix 1). 

Chromosomal partition: The analysis of chromosomal data resulted in 36 informa- 
tive characters, 23 constant characters and 40 non-informative characters. After ana- 



Sex chromosome systems in Alouatta 



55 




r 

I A. 



C. libidinosus 


2N 


54. XX/XY 9 (? 




L. lagothricha 


2N 


= 62, XX/XY 9 ( ? 




A. pigra 


2N 


= 58,X,X,X ! X ! /XiX2Y,Y 2 9<J 




A. pattiata 


2N 


= 54, X,X,X 2 X 2 9 and 53, X,X 2 


Y<? 


A. belzebul 


2N 


= 50,X 1 X,X 2 X 2 9and 2N=49, 


X,X 2 Y g 


A. caraya 


2N 


= 52, X^.X^XiXjYJj 9 (J 




A. macconnelli 


2N 


= 47 to 49, X.X.X^/X.XjY.Y;, 


9<?(m) 


A. sara 


2N 


= 48 to 51 , x,x,x 2 xjx,x 2 y,t 2 


9 <? (m) 


A. seniculus arctoidea 


2N 


= 44 to 47, X,X,X 2 X 2 /X,X 2 Y,Y 2 


9<J(m) 


A. g. clamitans 


2N 


= 46, X 1 X 1 X 2 X 2 X 3 X 3 9 and2N 


= 45, X, 


A. g. guariba 


2N 


= 50, X^.XjXjXjXj 9 and 2N 


= 49, X, 



Figure 4. 50% majority consensus tree obtained by "bootstrap" for the combined analysis. Next 
to the name of each species, the diploid number (2N) and sex chromosome system is described. 
(m)=microchromosomes. 



lyzing 704 trees, PAUP retained the two most parsimonious trees (Appendix 2: Figures 
Sa and Sb), both with a length of 87 (L = 87). The analysis using only the partition of 
chromosomal data did not resolve the node {{A. palliata, A. pigra), (A. caraya, A. bel- 
zebul), {{A. g. clamitans, A. g. guariba), (A. macconnelli (A. sara, A. s. arctoidea))), since 
it was established in a polytomy (Appendix 2: Figure Sc). 

Molecular partition: Heuristic analysis of cyt b gene sequences, made from a total 
of 800 characters, produced 109 informative characters, 551 constant characters and 
140 non-informative characters. After analyzing 916 trees, PAUP retained a single 
most parsimonious tree (Appendix 2: Figure Sd), with a length of L = 366. The analy- 
sis using only molecular data did not resolve the node {A. sara, A. macconnelli, A. s. 
arctoidea, A. caraya), which was established as a polytomy different from that described 
from chromosomal data. 

Combined analysis: The heuristic analysis of the combined data showed a total of 
899 characters, 145 of which were informative, 180 non-informative and 574 con- 
stant. After analyzing 684 trees, PAUP retained only one, with a length of L=460 
(Figure 4). This type of analysis allowed us to solve all the nodes, resulting in a fully 
resolved tree. 



Discussion 

We present the first phylogenetic study using a combined analysis of chromosomal and 
molecular characters in Ceboidea to contribute to the characterization of the specio- 
genic processes in howler monkeys. The homoplasy distribution is likely to be different 
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in each dataset because these are subject to different constraints. Therefore, when dif- 
ferent datasets are analyzed simultaneously, the signal common to all of them is more 
likely to overwhelm the homoplasy signal on the data (Kluge 1989). 

In primates, few studies have compared and taken into account more than one type 
of character. Bonvicino et al. (2001) superimposed chromosomal information on the 
phylogeny obtained from molecular characters. Villalobos et al. (2004) used numerical 
and metric values that describe the karyotype, such as diploid number (2N) and funda- 
mental number (FN), in a combined phylogenetic analysis with morphological charac- 
ters. However, these values (2N, FN, etc) can be identical simply by chance and, if in- 
terpreted in a phylogenetic context, may be spurious indicators of relatedness (Dobigny 
et al. 2004). Our encoding strategy (using the rearrangements as characters) is quite 
similar to that used for morphological data but in cytogenetics one can retrieve informa- 
tion on the mutational event itself, something that is clearly not available to morpholo- 
gists. As such, chromosomal mutations that accumulate along the tree are comparable to 
transitions, transversions, and insertions/deletions in molecular phylogenies (Dobigny 
et al. 2004). Our combined phylogeny evidences the accuracy of this encoding strategy. 

In all the above-mentioned contributions, the X 1 X 1 X,X 2 /X 1 X 2 Y 1 Y 2 sex chromo- 
some system was proposed as the ancestral condition for the genus. However, as dis- 
cussed by Solari and Rahn (2005), the X^X^/X^Y sex chromosome system is 
simpler and is present in other genera of Neotropical Primates, such as Aotus Illiger, 
1811, Callimico Miranda Ribeiro, 1912, and Cacajao Lesson, 1840 (Ma et al. 1976, 
Seuanez et al. 1989, Moura-Pensin et al. 2001). The X^X^/X^Y sex chromosome 
system as an ancestral state appears to be a more parsimonious hypothesis. Moreover, 
since Mesoamerican howlers {A. pigra and A. palliatd) were poorly karyologically char- 
acterized at the time, data on these howlers are missing in all previous contributions. 



Homeology analysis 

The karyotypes of A. pigra and A. palliata share more syntenic associations with those 
of A. caraya and A. belzebul than with those of the "A. seniculus group" {A. s. arctoidea, 
A. sara, A. macconnelli, denominated as such because they were once all subspecies 
of A. seniculus together with A. s. seniculus Linnaeus, 1766, and A. s. stramineus Hill, 
1962). This supports the basal grouping of the A. pigra- A. palliata Mesoamerican clade 
and the basal grouping of A. belzebul among South American howlers. 

The chromosomal comparisons showed that A. pigra and A. palliata conserved the 
HSA8/18 and HSA14/15 syntenies, considered ancestral for Platyrrhini (Stanyon et 
al. 2008), as well as the HSA10/16/10/16 syntenic association, ancestral for Atelidae 
(de Oliveira et al. 2002), but lost the HSA3/21 synteny, ancestral for mammals (Mul- 
ler et al. 2000, Muller 2006). 

According to our combined phylogeny, the HSA2/20 and HSA5/7/5/7 syntenic 
associations, previously considered as synapomorphies of the A. caraya-A. belzebul 
group (de Oliveira et al. 2002), would be homoplasies (parallelism). The HSA16/4 
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syntenyc association would be ancestral for the genus and might either be absent in A. 
pigra, A. guariba and the A. seniculus group or might constitute a parallelism among A. 
palliata, A. caraya and A. belzebul. 

Like the HSA3/21 synteny, the HSA3/15 syntenic association, involved in the 
sex chromosome systems in South American howlers, is not present in Mesoamerican 
ones. This syntenic association of human 3/15 chromosomal segments has been de- 
scribed in other Atelidae species such as Ateles geoffroyi Kuhl, 1 820 and Ateles belzebul 
hibridus Geoffroy, 1829, although not associated with the sex chromosome system 
(Morescalchi et al. 1997), but not observed in Lagothrix Geoffroy, 1812, and Brachyte- 
les Spix, 1823 (Stanyon et al. 2001, de Oliveira et al. 2005). This association has 
not been observed in other genera of Neotropical Primates such as Cebus libidinosus 
Spix, 1823, or Saimiri boliviensis boliviensis Geoffroy & Blainville, 1834 (Mudry et al. 
2001). Therefore, the HSA3/15 syntenic association either could be interpreted as the 
ancestral condition for the family Atelidae, where the association with multiple sex 
chromosomes would be an evolutionary novelty (apomorphy) in howlers and the loss 
of the association a apomorphy for the Lagothrix and Brachyteles group, or could have 
appeared independently in Alouatta, involved in the Y- autosome translocation, and in 
Ateles, not involved in the sex chromosome system (de Oliveira et al. 2005). However, 
our results suggest that the HSA3/15 syntenic association is not an ancestral condition 
for Alouatta, since the most basal species (see Figure 4) A. pigra and A. palliata (this 
contribution) and A. belzebul (Consigliere et al. 1996) do not possess this association. 



Possible origin of the multivalents 

Taking into consideration the data obtained, a hypothesis can be proposed regarding 
the origin of the sex chromosome systems in the genus. Within the family Atelidae, 
with the exception of Alouatta, all genera have an XX/XY sex chromosome system. 
Therefore, it can be considered that the Alouatta ancestor possessed a chromosomal 
sex determination XX/XY, prior to the biogeographic separation of Mesoamerican and 
South American groups (see below). After this separation, both groups independently 
acquired the multiple sex chromosome systems currently observed through independ- 
ent Y-autosome translocations. 

The sex chromosome system XjX^X^/XjX^Y would have arisen independently 
in the lineages of Meso and South American howlers by a Y-autosome translocation 
(Figure 5a). In males, two fissions, one in Yp tcr and another in q prox of the autosomal 
pair involved (Aq ), followed by translocation of Aq to Yq-p , formed the new 

f iprox' ' iprox 1 f prox 

chromosome Y . The Yp segment is lost and the proximal region of the fissioned 
autosome either is lost or, in certain howler species, could have given rise to micro- 
chromosomes (e.g.: A. seniculus (Yunis et al. 1976, Lima and Seuanez 1991, Torres 
and Leibovici 2001), A. sara (Minezawa et al. 1985) and A. macconnelli (Lima et al. 
1990)). The homologous autosomal pair involved in the translocation is the one now 
denominated X,. In the case of South American howlers, the autosomal pair involved 
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Figure 5. a Possible origin for X X Xp^/X X,Y sex chromosome sysrem in rhe genus Alouatta. The 
ancestral X chromosome is shown in white, the ancestral Y chromosome in light gray and the autosomal 
pair (A) in dark gray. Two fissions occurs, one in Yp ter and another in q proi . of the autosomal pair involved 
(Aq prox ). The translocation of Yq-p to the Aq formed the new Yj chromosome and the homolog of the 
autosomal pair involved in the translocation is now denominated X_. The Yp acentric fragment is lost 
and the rest of the autosome (Ap and Aq pro;i ) could either be lost or remain as a microchromosome in some 
howlers b Possible origin for the X 1 X I X 2 X_/X 1 X 2 Y 1 Y 2 sex chromosome systems from a X ] X ] X,X 2 /X 1 X 2 Y 
system. The ancestral X is shown in white, the ancestral Y in light gray, the autosomal pair involved in the 
first translocation (A) in dark gray and the autosomal pair (A') involved in the formation of this new sex 
chromosome system in black. Simultaneous breaks in Xjp and A*q followed by the translocation of 
the rest of the A' q to X 2 p pro> . give origin to the new Xj chromosome. The X 2 p ter acentric fragment could be 
lost and the rest of the autosome (A') could either be lost or remain as a microchromosome in some howlers. 
The homolog to the autosomal chromosome in question is now identified as Y 2 c Simultaneous breaks in 
Y q and Y 2 q and a translocation between Y and Y 2 further explain the hybridization pattern observed in 
the sex chromosome systems of South American howlers. A de novo centromere arises in the remains of 
the old Yj (now YJ. The remains of the old Y 2 could either be lost or remain as a microchromosome in 
some howlers d Hybridization pattern in South American howlers. 
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would share homeology with HSA3, whereas in the Mesoamerican species it would 
share homeology with HSA7. 

From this XjXjXpC/X JCY sex chromosome system, an X^X^/X^Y ' Y ' sys- 
tem could have arisen from a new translocation (Figure 5b). Under this hypothesis, 
simultaneous breaks in X n p and q of another autosome (A'q ), followed by the 

z-t prox J-prox v iprox' J 

translocation of most of the A'q arm to X n , gave rise to the new X, chromosome. 

i It prox O 2 

The X n p acentric fragment is lost and the rest of the autosome (A') either is lost or 
could have remained as a microchromosome in some howler species (see above). The 
chromosome homologous to the autosome in question (A v ) became Y . In the case of 
South American howlers, the new autosomal pair involved in the sex chromosome sys- 
tem would share homeology with HSA15. A further translocation between Y and Y 2 
(Figure 5c) would explain the hybridization pattern of the segments with homeology 
to human chromosomes 3 and 1 5 observed in the sex chromosome systems X^X^X 2 X 2 / 
XjX 2 YY 2 in South American howlers (Figure 5d). 

On the other hand, in the Mesoamerican species, the X 1 X I X 2 X 2 /X ] X 2 Y 1 Y 2 sex 
chromosome system could have arisen either as described in Figure 5b (with the auto- 
somal pair involved sharing homeology with a human chromosome not yet identified 
by G-banding pattern) or by a fission in Y that would have given rise to two chromo- 
somes, the new Y (containing the segment corresponding to the ancestral Y chromo- 
some) and Y 2 (containing a portion of the autosomal pair with homeology to HSA7). 
This last hypothesis would require a centromeric activation in Y . 

However, considering the observation of the independent origin of the multiple sex 
chromosome systems in these two groups of howlers, the possibility of an independent 
origin of the X t X X^XJX X_Y and X X X^/X X^Y Y, sex chromosome systems within 
the Meso and South American groups cannot be ruled out until further studies. 

It can be considered that multiple sex chromosome systems would be an extremely 
rare phenomenon due to complication in meiosis. Extreme cases are platypus and 
echidna, with a large number of sex chromosomes (Bick and Jackson 1967, Renz et 
al. 2007). In primates, multiple sex chromosome systems are even more infrequent. 
Moreover, Alouatta would be the first case where an independent origin of multiple sex 
chromosome systems is described. In other taxa, such as Drosophila (Flores et al. 2008), 
Erythrinidae fishes (Cioffi et al. 2013) and mole-rats (Deuve et al. 2006), a few cases 
have been observed, but these descriptions are still scarce. 



Phylogeny of Alouatta 

The chromosomal homeologies and FISH analysis were used to construct a data ma- 
trix for the phylogenetic analysis. For comparison purposes, independent phylogenetic 
reconstructions were performed with each type of partition (Appendix 2: Figure Sa, b, 
c and d), along with the combined analysis of the two datasets (Figure 4). The chromo- 
some partition grouped A. caraya and A. belzebul as sister taxa, in agreement with that 
reported by de Oliveira et al. (2002), a relationship that was not observed in the other 
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two analyses, which grouped A. caraya with A. sara, A. s. arctoidea and A. macconnelli 
(although in the case of the molecular partition this relationship constituted a polyto- 
my). This last species arrangement was also proposed by Nascimento et al. (2005) and 
Cortes-Ortiz et al. (2003) using molecular characters. In our molecular data partition 
(Figure Sd), A. belzebul was grouped with the clade of A. guariba, in agreement with 
that reported by Bonvicino et al. (2001) and Cortes-Ortiz et al. (2003). The three types 
of analyses agreed to place A. sara, A. s. arctoidea and A. macconnelli into a single group, 
although the molecular data partition did not resolve the relationships between them, 
as they formed a polytomy. The grouping of all those species ("A. seniculus group") was 
observed in all phylogenetic studies performed so far (see above). Another coincidence 
was that the Mesoamerican species were placed as a separate clade from other South 
American species and, as expected, the two subspecies of A. guariba in one group. All 
sets of taxa analyzed in this new approach were solved without polytomies only with 
the combined analysis, demonstrating the usefulness of incorporating more than one 
source of data for a more accurate elucidation of the relationships among current taxa. 

The grouping of South American species as a separate group of the Mesoameri- 
can group coincides with previous phylogenetic analyses using only molecular charac- 
ters (Cortes-Ortiz et al. 2003, Ellsworth and Hoelzer 2006) and with the hypothesis 
of monophyletic origin of the Mesoamerican howlers previously proposed by Smith 
(1970). Smith's hypothesis holds that Mesoamerican howlers originated by an expan- 
sion of the geographic distribution of South American howlers after the formation of 
the Isthmus of Panama, estimated to be completed about 3 million years ago (Coates 
et al. 2004). However, other studies indicate that the rise of the isthmus was a process 
rather than an event (Knowlton and Weigt 1998), resulting in intermittent periods 
with connected and divided lands during the past 18 million years. Given this last 
fact, another hypothesis was postulated to explain the current geographic distribution 
of the species of the genus. Instead of a single colonization event, various founder 
events, either across the Isthmus of Panama during one of the periods in which the two 
Americas were connected or across islands in the Caribbean archipielago, could have 
occurred (Cortes-Ortiz et al. 2003, Ellsworth and Hoelzer 2006). Primate fossils have 
been found in Cuba and Jamaica, but the origin and relationships of these specimens 
with modern Platyrrhini are still under debate (Fleagle 1999, Gutierrez Calvache and 
Jaimez Salgado 2007, Rosenberger et al. 2009, Cooke et al. 201 1). 

Independently of the biogeographic scenario under consideration, it is clear that the 
evolutionary history of Mesoamerican howlers is different from that of South American 
howlers, an assertion that would be supported by the evidence provided by our new data. 

This contribution provides new useful information for the systematics of the genus 
Alouatta, while supporting the hypothesis of chromosomal evolution in primates as a 
speciogenic strategy. The combined analysis resolved the phylogenetic relationships 
between howler species of both American origins, as a first approach to the "Total 
Evidence" concept and towards clarifying the controversies related to the Taxonomy 
and Evolution of Ceboidea. 
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Appendix I 

Data matrix, (doi: 10.3897/CompCytogen.v8il.6716.appl) File format: 
Microsoft Word file (doc). 

Explanation note: Data matrix contains: Chromosomal data matrix, Molecular data 
matrix and Combined data matrix. 

Copyright notice: This dataset is made available under the Open Database License 
(http://opendatacommons.org/licenses/odbl/LO/). The Open Database License 
(ODbL) is a license agreement intended to allow users to freely share, modify, and use 
this Dataset while maintaining this same freedom for others, provided that the original 
source and author(s) are credited. 
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Appendix 2 

Supplementary Figure S. (doi: 10.3897/CompCytogen.v8il.67l6.app2) File 
format: Microsoft Word file (doc). 

Explanation note: Figure S: a) and b) Most parsimonious trees obtained for the chro- 
mosomal partition c) 50% majority consensus tree obtained by "bootstrap" d) 50% 
majority consensus tree obtained by "bootstrap" for the molecular partition. 

Copyright notice: This dataset is made available under the Open Database License 
(http://opendatacommons.org/licenses/odbl/LO/). The Open Database License 
(ODbL) is a license agreement intended to allow users to freely share, modify, and use 
this Dataset while maintaining this same freedom for others, provided that the original 
source and author(s) are credited. 
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